manuscript no. isolat 


February 2, 2008 


(DOI: will be inserted by hand later) 





o 
o 

> 

o 



> 
cn 
cn 
in 
a^ 
o 

o 
6 



A Solution to the Isolatitude, Equi-area, Hierarchical 

Pixel-Coordinate System 

Boudewijn F. Roukema, Bartosz Lew 

Torun Centre for Astronomy, N. Copernicus University, ul. Gagarina 11, PL-87-100 Torun, Poland 
Le 2 fevrier 2008 

Abstract. Gorski et al. (1999b) have earlier presented the outline of a pixelisation-to-spherical-coordinate trans- 
formation scheme which simultaneously satisfies three properties which are especially useful for rapid analyses of 
maps on a sphere: (i) equal spacing of pixels along lines of constant latitude, (ii) equal pixel 'areas' (solid angles) 
and (iii) hierarchical scaling with increasing numbers of pixels. Their outline is based on the division of the sphere 
into twelve regions covering equal solid angles, which are hierarchically subdivided in a way compatible with 
these three criteria. In this paper, a complete derivation of this scheme is presented, including, in particular, (1) 
the angle 6* defining the limit between polar and equatorial regions, and (2) the transformations from the unit 
interval [0, 1] A [0, 1] to spherical coordinates in a polar region. 
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1. Introduction 

The most commonly required analyses of microwave back- 
ground data require spherical harmonic analyses of tem- 
perature fluctuation maps. Given increasingly large num- 
bers of pixels, Gorski, Hivon, & Wandelt (1999a); Gorski 
et al. (1999b) pointed out that an ideal pixelisation-to- 
spherical-coordinate transformation scheme would satisfy 
the following requirements: 

(i) an isolatitude requirement, i.e. pixels should be equally 
spaced in latitude, 

(ii) pixel solid angles (hereafter, 'areas') should be equal, 
and 

(iii) the system should scale hierarchically with increasing 
numbers of pixels, preferably according to some power 
of two, in order to maximally exploit the binary nature 
of present computer technology. 

Gorski et al. (1999b) present the outline of a system 
satisfying these requirements. The idea of dividing the 
sphere into 12 regions — four equatorial regions and eight 
polar regions — which are mapped to or from the unit 
interval [0, 1] A [0, 1], is presented there. 

However, the full solution is not derived in Gorski et al. 
(1999b). 

In this paper, firstly, in Sect. 12.11 the angle 6* defin- 
ing the limit between polar and equatorial regions of this 
scheme is derived. 




4 half polar regions 
1/6 of full sky =2Ji/3ster 

4 half polar regions 
+ 4 half equatorial regions 
= 1/3 of full sky 

4 half polar regions 
+ 4 half equatorial regions 



4 half polar regions 



Fig. 1. Schematic diagram of the 12 regions, or equivalently, 
24 half regions, showing how the isolatitude requirement and 
the equal pixel area requirement together require that the 
solid angle from the line of constant latitude 6* to the pole 
is 2-k/3 ster. 



As is described in Gorski et al. (1999b), the binary bits 
of fractions with the square unit interval [0, 1] A [0, 1] are 
used for pixel ordering within a region. For completeness, 
we remind the reader of this ordering in Sect. 12.21 and 
express it algebraically. 

In order to have the complete transformation from 
pixel numbers to spherical angles, the transformation be- 
tween the fractional position {fi, fj) in the square unit 
interval and a pair of spherical angles (0, 9) is needed. In 
12.3.11 this transformation is derived for equatorial 
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for polar regions. 
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Fig. 2. Ordering of the 12 regions, shown in a projection where the fuU sky, 4n ster, is shown as 12 square unit intervals. Each 
pole is represented as a set of four distinct points in this projection, indicated by arrows. The derivation of the straight line (in 
this projection) dividing the "polar" regions into "polar" and "equatorial" parts, labelled sin6' = 2/3, is derived in Eq. O. 



The convention chosen here for spherical coordinates 
is to have a longitude coordinate 

0e[O,2^] (1) 

and a latitude coordinate 

0e [-V2,7r/2]. (2) 

For clarity of presentation, it is also useful to define 
(jjT, = (p/n (3) 
and 

s = sine. (4) 
2. Calculations 

2.1. The boundary between equatorial and polar 
regions: 9* 

Fig.n(cf. fig. 2 of Gorski et al. (1999b)) and Fig.Hshow 
the overall geometry of the 12 regions, including an order- 
ing, though any other ordering could easily be chosen. 

In order for the 12 regions to have equal pixel areas, 
the coordinates of the four vertices along the equator, of 
the equatorial regions, must be spaced by an interval of 
A(/) = 7r/2, i.e. they are 



The two polar vertices are 

0^ undefined, 6* = -7r/2,+7r/2. (6) 

Finding the coordinates of the remaining eight vertices, 
which separate equatorial and polar regions, requires an 
integration. 

Define 9* to be the angle from the equator to the curve 
of constant latitude separating the equatorial regions from 
the "polar" half of the polar regions, as shown in Fig. ^ 
(The "polar" regions are, in fact, divided between a "po- 
lar" half which is polar and an "equatorial" half which 
uses a transformation similar to that of the "equatorial" 
regions. This is explained below in Sect. 12. ^T^ ) 

Since the system is an isolatitude, equal pixel area sys- 
tem, the entire solid angle from a curve of constant lati- 
tude 9* to the pole must be half of the total solid angle 
represented by four regions, i.e. half of one third of the 
full sky, i.e. one sixth of Att steradians, i.e. 27r/3 ster. This 
is illustrated in Fig. ^ 

Hence, 

/ d0 / cos 6* de = 27r/3, 

^ 27r[sin6']g/^ = 27r/3 
27r(l -sinr) = 27r/3 



0, = 1/4, 3/4, 5/4, 7/4,^ = 0. 



(5) 



=^sinr = 2/3 



(7) 
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2.2. Hierarchical ordering of pixels 



2 sine*=2/3 
1 equator 

sine*=-2/3 



Fig. 3. Equatorial region: unit interval [0, 1] A [0, 1] 
parametrised by fractions fi, fj and linearly transformed to 
X = fi — fj and y = fi + fj. This can be any of the regions 5, 
6, 7, 8; without loss of generality, the discussion in Sect. 12.3.11 
is presented for region 5. 




2 pole 

1 sine*=2/3 

equator 



Fig. 4. Polar region: unit interval [0, 1] A [0, 1] parametrised 
by fractions fi , fj and linearly transformed to x = fi — fj and 
y = fi + fj. This can be any of the regions 1, 2, 3, 4, 9, 10, 
11, 12; without loss of generality, the discussion in Sect. 12.331 
is presented for region 1. 



^0* = sin-i 2/3 « 41.8° 



(8) 



Hence, in Fig.|21 the curve of constant latitude dividing 
each polar region into two equal halves is labelled sin 9 = 
±2/3. 

Note that 6* is not 7r/4, even though a casual glance 
at Fig. ^a-iid casual intuition might suggest this. 

Hence, the eight vertices separating equatorial and po- 
lar regions are 



= 0, 1/2, 1, 3/2, d = r), 
0, 1/2, 1, 3/2, 6 = -r). 



(9) 



Requirement (iii), the hierarchical ordering, is a binary or- 
dering within the (square) unit interval, and an ordering 
of the twelve unit intervals. This ordering is referred to as 
'nested' by Gorski et al. (1999b), and is the ordering used 
in the data files of the first year WMAP data (Bennett 
et al. 2003) , where and 9 are galactic longitude and lat- 
itude respectively. It is shown for two different resolution 
levels in the lower half of fig. 3 of Gorski et al. (1999b). 

2.2.1. Binary ordering of pixels within the square unit 
interval 

Figs H show the definit ions of fi , fj as fractions in the 
unit interval in an equatorial and a polar region respec- 
tively. 

Let us write fi, fj as binary fractions, rounded to the 
nearest pixel centre, to some precision n > 1, correspond- 
ing to linear pixel sizes of approximately, but not exactly, 

it/2 rad -1 f 41 7r/3 ster 
' ^„ — , or pixel areas oi exactly ' ^„ : 

0<fi= Q.b2n-l-b3bil, fj = 0.b2n-b4b2l < 1 (10) 

where bi G {0, 1}. The last binary place is occupied by "1" 
in order to have a pixel centre. Then the bits in fi and 
fj provide alternating bits for integer ordering of pixels 
within the square unit interval. 

This mapping is sufficient for both equatorial and polar 
regions. 

This is shown in fig. 1 of Gorski et al. (1999b). 
For example, with just one binary digit, a region is 
divided into four pixels. If we have binary representations 



/, =0.6il,/, =0.621, 



(11) 



then the pixel centred at (fi, fj) is the (626i)2-th pixel 
within the region. 

With two binary digits, a region is divided into sixteen 
pixels. If we have binary representations 



/, = o.Mii,/j =0.64621, 



(12) 



then the pixel centred at {fi, fj) is the (6463626i)2-th pixel 
within the region. 

More generally, with n > 1 binary digits, the pixel at 

ihJj) 

= 0.62„-1...636il, /, = 0.62n...6462l, (13) 

is the (62n---626i)2-th pixel (again in binary) in the region. 

2.2.2. Ordering of the twelve square unit intervals 

The ordering of the twelve regions shown in Fig. [21 is that 
used in the first year WMAP data files. This combines 
with the ordering within a region. 

That is, with n > 1 binary digits, the pixel centred at 
ifiifj) in the k-th region, with binary representations 



fi = 0.62„_1...636il, /, = 0.62n---6462l, 



(14) 



Roukema et al.: Isolatitude, Equi-area, Hierarchical Pixelisation 



Table 1. Vertex conversions for an equatorial region. 
Equations (|16|l . ^ and define x,y,(f>Tr and s. See eqns 
JSJ, © and lO for the values of (/i^ and s. 



label 


X 


y 




s 


A 











-2/3 


B 


1 


1 


1/4 





C 





2 





2/3 



(17) 



is pixel number 

[(fc-l)22"]^^ + [62„...626i]2 (15) 

where the first term is shown in decimal notation and the 
second in binary notation. 

2.3. Transformations from the unit interval to 4>, 9 

For either an equatorial or polar region, shown in Figs|31 
01 respectively, it is useful to make the change of variables 



y 



fi fj 



(16) 



where /i, fj are the fractions in the unit interval. 
The requirements of 

(i) isolatitude pixel positions, and 

(ii) equal pixel areas 

then require transformations from {fi, fj) to {<j),0), or 
equivalently from [x, y) to (0, 0), or even more simply from 
[x,y) to ((/)7r,s), as defined in eqs ©, l@J. 

Use of equal intervals in s = sin 6 in latitude enables 
straightforward linear mappings of the complete equato- 
rial band (the regions labelled here as "equatorial", plus 
half of each of the regions labelled here as "polar" ) . This 
can be thought of as stretching the latitude coordinate 
in order to compensate for shorter lengths of curves of 
constant latitude. 

However, this cannot function close to the poles, hence, 
a distinct mapping is required for the polar regions. 

2.3.1. Equatorial regions: transformations from unit 
interval to 6, 9 



Since 
d(sin6 



cos6i de. 



(18) 



the mapping in any of the four equatorial regions is pro- 
vided by a linear map from x to (p and from y to sm9. 

Without loss of generality (there are four equatorial 
regions in total), it is sufficient to consider the region la- 
belled "5" in Figs ^ El to be the region which has </) = 
at its centre. 

Then, the definitions ((TB|l . (jSJ and Q together with 
the values in eqns lO), lO and ^ give the vertex conver- 
sions shown in Tabled 

The requirements of 



Table 2. Vertex conversions (B to E) and a boundary 
condition (F) for a polar region. Equations IjlBp . and 
Q define x,y,(j)Tr and s. See eqns © and for the 
values of (^^r and s. 



label 



X y 



B 
C 
D 
E 

F 




-1 1 
-fl 1 

2 



< lim 



1/4 
2/3 
1/2 2/3 



(20) 



r{-x = y - 2) = <f,^{x = y - 2) + 1/2 



(i) isolatitude pixel positions, and 

(ii) equal pixel areas 

together with the vertex conversions in Tabled yield the 
linear transformations 



0^ = x/A 



(19) 



2.3.2. Polar regions: transformations from unit interval 

to 4),e 

Without loss of generality, the polar region "1" , as shown 
in Figsnm adjacent to the equatorial region labelled "5", 
is considered. 

The vertex conversions for this polar region are those 
in Table El 

Vertex conversion Eq. H2U|) .F is a requirement of conti- 
nuity of between adjacent polar regions: the two bound- 
aries of the upper half of the region are the lines x = y — 2 
and —x = y — 2 (respectively upper left and upper right 
boundaries of the unit interval in Fig.QJ and for any fixed 
value of y, must cover the angular interval [0, 7r/2], i.e. 
the 0^ value at constant y must increase by A^i^ = 1/2 
between these two limits. Hence, Eq. H2Uf) .F. 

The region is illustrated in Fig. 01 

The lower half of the 'polar' region, where y < 1, is 
really in the equatorial part of the sphere, so a similar 
linear solution applies as above: 

l + x 



2y 
3 ' 



(21) 



providing continuity and satisfying the isolatitude and 
equi-pixel requirements and the vertex conversions (|20() .B. 
(113). C and (Ill.D in TableEl 

By continuity, I^HI .C and (jSOJ-E) must also be satisfied 
for limiting cases in the upper half of the region as y ^ 1+. 

The isolatitude and equi-pixel requirements imply that 
for y > 1, 4>Tr is linear in x, and in order to satisfy (|20fl .C 
and (EOJ-I^, it follows that 

9^ = (22) 



Roukema et al.: Isolatitude, Equi- 

for some function / which does not depend on x and which 
satisfies 



fiy = 1) = 1 

in order to be continuous with the region y < 1. 
Boundary condition Eq. ^^.F (Table EJ is 

(P^i-x = y-2) = (l>^{x^y-2) + l/2 
_ l + /(y)[-(y-2)] _l + f{y){y-2) 



4 4 
/(y)h(y-2)] = /(y)(y-2) + 2 

- 2^y 



1/2 



(23) 

(24) 

(25) 
(26) 
(27) 



Because boundary condition F implicitly includes con- 
ditions C and D when y = f, the required constraint in 
Eq. ^ is satisfied, i.e. /(I) = 1. 

Thus, 



2-y + x 
4(2 - y) ■ 



(28) 



To obtain s, note that equality of pixel areas for y > 1 
and y < 1 requires that 



(29) 



ds 






ds 






dy 


y>l dx 




dy 


y<l dx 


y<l 



ds 
dy 



Equations HU then give 

1 _ 2 f 
4(2^ "34 

2 



ds_ 
dy 



^(2-y). 



(30) 
(31) 



y>l 



The isolatitude condition implies that 

ds ^ 
dx ' 

so integration gives 

s = constant — -{2 — y)^, 
3 



(32) 



(33) 



and vertex conversion D provides the constant, so that 



(34) 



Taking the limit of as {x,y) — > (0,2), within the 
range y — 2 < x < 2 — y, clearly provides vertex conversion 
E. 
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3. Conclusion 

The solution to the regionalisation of the sphere recom- 
mended by Gorski et al. (1999b), with the three require- 
ments 

(i) isolatitude, i.e. equal spacing of pixels in latitude, 
along lines of constant latitude, 

(ii) equal pixel areas, and 

(iii) hierarchical ordering of pixels within square unit in- 
tervals which map to equal solid angles, 

has been completed by deriving, in particular, the angle 9* 
separating polar and equatorial regions, and the full trans- 
formations from the unit interval [0, 1] A [0, 1] to spherical 
coordinates in a polar region. 

Together with the binary ordering presented by Gorski 
et al. (1999b), the full solution is as follows. 

Binary, bitwise mappings relate fractional positions 
ifii fj) within the fc-th square unit interval to a pixel num- 
ber. With n > 1 binary digits, the pixel centred at {fi, fj) 
with binary representations 



fi = 0.62«-l---&3&ll, /l = 0.62„...64&2l, 



(35) 



where bi £ {0,1}, in the fc-th region shown in Fig. El is 
pixel number 



[(fc-l)22"]^^-f [62„...526l]2, 



(36) 



where the first term is shown in decimal notation and the 
second in binary notation, as mentioned above in Eq. (|15|l . 

The vertices along the equator of the equatorial regions 
are (from Eq. ©) 



0^ = 1/4, 3/4, 5/4, 7/4,0 = 0; 



(37) 



the vertices separating equatorial and polar regions are 
(from Eq. ®) 

(</),= 0, 1/2, 1, 3/2, = 0*), 

((/), =0, 1/2, 1, 3/2, = -r); (38) 

and the two polar vertices are (from Eq. lO) 

(f)^ undefined, = -7r/2, +7r/2. (39) 

Mappings from {x = fi - fj,y = fi + fj) to {(j)^ = 
(/)/7r, s = sin0), are found, without loss of generality, for 
one equatorial region and one polar region: 



j-^^** ^-^ \ equatorial region 



_ l+x 

I 4 

3 



/ _ 2-y+x 

. = l-i(2-y)2 



polar region, y < 1 

polar region, y > 1 (40) 



summarised here for convenience from equations (|19|l . 
(Pljl. (PI) and 
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